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Gels of semi-flexible polymers, network glasses made of low valence elements, softly compressed ellipses and dense 
suspensions under flow are examples of floppy materials. These systems present collective motions with almost no 
restoring force. We study theoretically and numerically the frequency-dependence of the response of these materials, 
and the length scales that characterize their elasticity. We show that isotropic floppy elastic networks present a phonon 
gap for frequencies smaller than a frequency uj* governed by coordination, and that the elastic response is localized 
on a length scale l c ~ 1/y/uJ* that diverges as the phonon gap vanishes (with a logarithmic correction in the two 
dimensional case). l c also characterizes velocity correlations under shear, whereas another length scale I* ~ 1/uj* 
characterizes the effect of pinning boundaries on elasticity. We discuss the implications of our findings for suspensions 
flows, and the correspondence between floppy materials and amorphous solids near unjamming, where l c and I* have 
also been identified but where their roles are not fully understood. 



1 Introduction 

In 1864 Maxwell^ showed that in order to be mechan- 
ically stable, the average coordination z of structures 
made of points connected by rigid bars must be larger 
than a threshold value z c . Most common solids satisfy 
a microscopic version of this constraint.'^ If the con- 
straint is violated, collective modes with no restoring 
forces (floppy modes) exist. However, in some cases these 
modes arc stabilized by weak interactions: the bending 
energy of semi-flexible polymers confers a finite elastic- 
ity to gels,^ as do van der Waals interactions in weakly- 
coordinated covalent glasses.^ Floppy modes can also be 
stabilized by anharmonicities, as is the case for gently 
compressed packings of ellipses.^ Fluids, on the other 
hand, can display modes that are strictly floppy. For 
example, in granular systems or suspensions of hard par- 
ticles, large and sometimes percolating clusters of con- 
nected particles can be formed/^ an d motion within 
these clusters can only occur along floppy modes where 
no particles overlap. As the density increases toward jam- 
ming, the viscosity^^ and the length scale^HSl charac- 
terizing the correlation of the dynamics diverge, up to 
the point where floppy modes disappear and the dynam- 
ics stops. 

In ellipses QlilASJ models of covalent networks,^ gels of 
semi-flexible polymers^ and suspension flows^ it has 
been observed that the density of vibrational modes D(u>) 
displays a gap at low frequency (in addition to the floppy 
modes present at zero-frequency) . In suspension flows the 
amplitude of the gap was shown to affect the divergence of 
the viscosity near jamming.^ An early work by Garboczi 
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and Thorpe®! supported the idea that a gap of modes 
exist in floppy materials. However, the dependence of the 
gap on the microscopic structure, and its consequences 
on the material properties and the different length scales 
characterizing the elastic response are not understood. 

In the present paper, we study elastic networks that 
are strongly disordered, but where spatial fluctuations of 
coordination are weak. The model networks we use are 
presented and analyzed numerically in section [2| In sec- 
tion [3J we show, using both mean-field methods^H^H and 
numerical simulations, that such isotropic floppy elastic 
networks present a vibrational gap below a frequency 
to* ~ z c — z = Sz. We show also that the absence of 
low frequency modes suppress the elastic propagation of 
the response to a local strain, which is thus localized 
and shown to decay exponentially as exp(— r/l c ), where 
l c ~ l/y/5z. In sectionUJwe go beyond mean field and es- 
timate the fluctuations around the mean response, which 
dominate the amplitude of the response and decays with 
the same length scale ~ l c . We also show that although l c 
can be considered to be the localization length of floppy 
modes, these modes can be strictly localized (have a com- 
pact support) only on the much larger scale I* ~ 1/Sz. 
In section [5] we study the floppy networks stabilized by 
weak interactions, which confer a finite but small restor- 
ing force to the material. We find that localization is 
lost even when weak interactions have a vanishingly small 
amplitude, but l c still characterizes the response near the 
applied strain. In section [6] we compare our results with 
recent observations in the affine solvent model of suspen- 
sion flows, where the spectrum associated to the evolution 
operator shares both commonalities and striking differ- 
ences with the density of states of the isotropic floppy 
networks considered here. Finally, our results raise ques- 



tions associated with the respective role of l c and I* both 
in floppy materials and in jammed packings, which are 
discussed in the last section. 

2 Model description and simulation 

We consider a network of N point particles of mass m, 
connected by N c harmonic un-stretched springs in spatial 
dimension d. The elastic energy following a deformation 
vector field {<5R l }i=i...Ar is: 

5E = E ^[(^ " SKJ ) ■ n ^ 2 + o(5B. 2 ), (1) 

(ij) 

where the index (ij) runs over the N c springs, whose 
stiffnesses are fcjj and whose directions are along the unit 
vectors n^-. Floppy modes correspond to displacements 
for which SE = 0. Since the energy is the sum of N c pos- 
itive definite terms, it vanishes only if all these terms are 
zero. Such modes always exist if the number of degrees 
of freedom Nd is larger than the number of constraints 
N c , or equivalently if z = 2N C /N < z c = 2d. 

The force field |F) generated by displacements can be 
obtained by taking the derivative of Eq. ([!]). One obtains 
|F) = M\SR), where: 

M = ^hjxiij ® mj(\i) - \j))({i\ - 01) (2) 

(ij) 

is the stiffness matrix. We have used the bra-ket notation 
for which (i|5R) = 5Hi. The normal modes are given by 
the eigenvectors of M. , and their frequencies are given by 
the square root of their associated eigenvalues. Herein, 
we will investigate the spectrum of M. and the spatial 
properties of its eigenvectors when z < z c . 

As a model system we consider isotropic disordered 
elastic networks, generated by following the method of 
RefJ^Sl we prepare amorphous packings of compressed 
soft elastic particles with a coordination significantly 
larger than z c . The centers of the particles are the 
nodes of our networks, and the contacts between par- 
ticles are replaced by un-stretched springs of stiffness 
k. Springs are then removed until the desired coordina- 
tion z is reached. Removal takes place randomly from 
the set of most connected pairs of nodes, leadings to 
isotropic networks with low heterogeneity in density and 
coordination. Such networks are appropriate models of 
amorphous solids for which large spatial heterogeneities 
are not e nergetically favorableJ 23 * 25 ! In rigidity percola- 
^j on 120J22] S p rm g removal is completely random, resulting 
in networks with large fluctuations that affect the clastic 
properties. 



For these networks, we diagonalize M. numerically to 
compute the density of states D(lo) for various z < z c . 
As appears in the inset of Figjl] we find that a gap in 
the vibrational spectrum appears below some frequency 
uj* that decreases as z — > z c . If spatial fluctuations of 
coordination were large, one would expect this gap to be 
filled-up in the thermodynamic limit due to Griffiths-like 
singularities, but we do not observed this effect in the 
network considered here. The role of this gap in elas- 
ticity can be analyzed by considering the response to a 
local strain. We change the rest length of one spring and 
let the system relax to zero energy under over-damped 
dynamics. Such a response can also be obtained in the 
absence of damping, by imposing a local oscillatory strain 
at a vanishing small frequency. As shown in Fig[2j the 
elastic information does not propagate: the response is 
localized on some length scale that appears to diverge as 
the gap vanishes, i.e as z — > z c . The response appears to 
be very heterogeneous. In what follows we will investi- 
gate its mean and its fluctuations. 




Fig. 1: Rescaled density of states vs. rescaled frequency 
for z £ [3.2, 3.95] using N = 10000 nodes in two dimensions. 
The continuous line corresponds to the theoretical prediction 
of Eq.(|9|. Inset: Non-rescaled density of states D(w) vs u. 
Floppy modes lead to a delta function at cj — and are not 
presented. 



3 Effective medium theory (EMT) 
3.1 General formalism 

We use EMT, also known as coherent potential approxi- 
mation, ^ 1 * 22 * 26 ! to investigate the behavior of floppy net- 
works as the one observed in the previous section. EMT is 



a mean field approximation and neglects spatial fluctua- 
tions in the coordination. Therefore, it does not properly 
describe rigidity percolation.^ However, it was shown to 
describe materials for which spatial fluctuations of coor- 
dination are small.^ EMT attempts to describe disor- 
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Fig. 2: Displacement field caused by a local strain in a floppy 
network. An over-damped relaxation was used following the 
elongation of one spring for (a) Sz — 0.2 and (b) Sz = 0.05. 



dered materials as ordered materials with a frequency- 
dependent effective stiffness fc ff. We apply this tech- 
nique to isotropic lattices of coordination Zi n >> z c , 
where bonds are then randomly removed with a prob- 
ability (1 — p), so that the stiffness coefficients fey take 
the values or fc with a probability (1 — p) or p respec- 
tively, and the final coordination is z — pzi n . The Green's 
function G(ui) of the disordered floppy system is defined 
as [— moj 2 + M]G(uj) = —1. Go is the Green's func- 
tion of the effective medium, corresponding to the initial 
ordered lattice with an undetermined effective stiffness 
coefficient k e g. Standard calculations lead to the Dyson 
relation G = G + GoTGo,^ where the operator T can 
be expressed as an infinite series in terms of increasing 
numbers of interacting contacts: 



E 



■ W. 



G T 



The transfer matrix is found to be: 

- \j)){hj - k e «)({i\ 



■ (ij) 



(kl) 



01) 



l-(%-fe eff )ny((i|-01)G (|i)-|i)K 



(3) 

We seek an effective stiffness fc e ff that captures the av- 
erage behavior of the system, i.e. (G) = Go, where the 
average is taken over the disorder on the stiffness coeffi- 
cients fej,-. This condition leads to (T) = 0. In the EMT 



this constraint is approximated by (T/y\) = 0. Using 
standard identities for the Green's function on isotropic 
lattices,^ one can express this condition as:^ 



mcj 2 tr[G (V = 0,w)] = 



Sz + (zjn - 2d)k cS 
2(1 - fc off ) : 



(4) 



where tr[»] stands for the trace, and fc e ff = fceff /fc- 

Since we are interested in low frequencies ui « y^k/m, 
we approximate Go by its continuum limit and use a 
Debye cut-off qo- We introduce the bulk modulus K and 
the shear modulus fi of the ordered lattice without an 
effective stiffness, and use the approximation: 



G (r, w ) = ]T / 
„ Jo 



d d q 



iqr 



2ir d moj2 



c a fceff <7 ; 



■q«®q a - (5) 



The sum is taken over all polarization vectors q Q , where 
c p = jr{K + |yu) for pressure waves and c s — j^/i for 
shear waves. 



3.2 Scaling analysis near z — z c 

We now perform a scaling analysis of Eqs.(4j5) as z ~> 
z c from below. Eq.([5]) implies that tr[Go(0,cj)] = 

-k S! f ( -"/Tfc ) ) wriere the function / is independent of 
z and u). As u) — > 0, the elastic moduli of floppy 
systems vanishes and fc G g — > 0. Eq.Q then leads to 

mu 2 _ _ i 

II ferff 



lim,,,^n m ,f . = e, where e satisfies the equation 



e/(e) - Sz/2. 



(6) 



For spatial dimensions d > 3, Eq.([5| implies that /(0) is 
a positive constant, therefore in the regime 8z << 1 one 
finds 

e«W(0)). 

Using Eq.Q together with the assumption that 
| mw 2 /fceff | << 1 and |fc c ff| << 1 (which can be shown 
to be true a posteriori in the limit of Sz << 1 and 
u> « Wkjm) one finds: 



fceff 



-Sz + V- 8 /(°)(^m - 2d)muj 2 + Sz 2 
2(z in - 2d) 



(7) 



For d = 2, Eq.([5]) leads to a logarithmic divergence in 
the small e limit: 



V 

'2N 



(c; 1 + c ; 1 )iog(e) + r, 



(8) 



where T is a constant that depends on the elastic moduli 
and the Debye cut-off. The weak divergence of / does not 
modify the perturbation analysis performed for d > 3, 



except that /(O) must now be replaced by /(e) in Eq.Q. 
Only /(e) remains undetermined. 

The asymptotic value for e in the limit 8z — > can 
be obtained from Eq.M, which leads to e ~ iogf/ z ) an< ^ 
therefore /(e) ~ — log(Tz). This behavior is only valid for 
values of | log(<5z)| 3> 1, a difficult limit to observe empir- 
ically. Therefore, to compare our theoretical predictions 
with two dimensional numerical observations we compute 
e(Sz) by solving Eq.Q numerically using Eq.Q. We find 
that the relation e(8z) depends on the initial network via 
only one parameter A(c p , c s , r)F] To compare the theory 
and numerical simulations, we use A(c p , c s , T) as a fitting 
parameter. 

The mechanical response of floppy systems is there- 
fore fully determined by the Green's function in Eq.([5]) 
and the effective stiffness given by Eq.Q. We start by 
computing the density of states D(ui) using the relation 
D(lu) = ^Im[tr[G o (0,w)]]. For < u < lo* , where 
lo* = Sz =, we find that D(u) = 0. For 

y/Sf(e)( Zi? -2d)m' V ' 

lo > to* we obtain: 



a) 



b) 



J • u 


Ji ii: <y~~ 






//■//•/,, 
//■//■/ 


w a 

, .vww 

V \ \ A \ 




I 7 / ■ I ! I \. \ \\ \ * 



Fig. 3: Mean response following the elongation of a spring 
placed at the center, along the vertical axis, a) Theoretical 



prediction given by Eq. ( 10 \ using c s and c p of an hexago- 
nal lattice. Arrows show the displacement field and contin- 
uous lines are stream lines of the vector field, b) Average 
response over 6000 independent realizations of networks with 
N = 40000 nodes and 8z = 0.05. 



D(cj) 



2/(e)(z m -2d)m 

^2 



'1- 



(9) 



Thus the size of the gap scales linearly with Sz for d > 3, 
where /(e) /(0), whereas for d = 2, a logarithmic 
correction exists. According to Eq.Q, rescaling the fre- 
quency by lo* and D(u>) by f(e) 1 ^ should collapse the 
low-frequency part of the spectrum. FigjT] shows that 
the quality of the collapse is very good. For all consid- 
ered coordinations, we used the same fitting parameter 
A(c p ,c s ,r) = 1.3, which is fixed by this measurement. 

3.3 Linear response to a local strain 

The response to a local force G (w, r^) diverges as ui — > 
due to the presence of floppy modes. One must rather 
consider the response field due to an imposed displace- 
ment. Changing the rest length of a spring placed be- 
tween nodes i and k at a fixed frequency loq corresponds 



iujat 



to imposing a displacement (<5R^ — <5Rfc) • n = e 
n being the unit vector along the connecting bond, 
force required to impose such a displacement is 



with 
The 



\F) = 



2n • (G (wo) 0) - G a (uj , r ki )) ■ n 



(\k)-\i)). 



In the small frequency regime, the last expression is found 
to be \F) ~ k cS e Wot n(\i) - \k)). The magnitude of the 
force vanishes as loq — > 0, which is consistent with the 
existence of floppy modes. 



* The parameter A(c p , c s , T) 



log[^( C - 1 +c7 1 )] + 



The response at the zero frequency limit in two dimen- 
sions, at distances r = |r| larger than the typical springs 
length (i.e. rqo S> 1) is given by 



(<5R(r)) = g(r - r ; ) - g(r - r k ) 



(10) 



where 



g( r ) = 



V x (V x K (rJ -)n) -V(V ■ K {rJ -)n) 



and K is the modified Bessel function of the second kind, 
that behaves exponentially at long distances. Thus, the 
mean response to a local strain decays exponentially in 
floppy systems with a characteristic length l c ~ l/\/e- 
The mean perturbation induced by changing the rest 
length of a spring has a quadrupolar symmetry shown in 
Fig[3^. After averaging over 6000 realizations, we obtain 
good agreement with our theoretical prediction, shown in 
in Fig. gjb. 

The asymptotic solution in any dimension has an ex- 
ponential decay. Indeed, taking the angular average one 
gets 



6R m (r) = (6R(r)) 



,1/2 



c (d-l)/2 -r/l c 



{ich 



\(d- 



-l)/2 j 
(11) 

where the over-line stands for angular average. In two 
dimensional networks, rescaling the distance by l c and 
the amplitude by ^/e collapses the response to a local 
strain at different coordination values into a single curve, 
as shown in Fig|4^-b. 
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Fig. 4: a) Mean response <5il m vs. the distance r from the 
imposed strain with N = 40000 nodes, b) Rescaled Average 
displacement of the mean response 8R m vs. the rescaled dis- 
tance r/l c , using A(c p , c s ,F) = 1.3 as extracted from Fig.Q. 
c) The fluctuations around the mean response are character- 



ized by SRt, defined in Eq.( 16 1. SR t is plotted vs. the distance 
from a local strain r. d) 5R t vs. the rescaled distance r/l c , 
using A(c p ,c s ,F) = 1.3. 



The zero frequency limit can be extended to finite fre- 
quencies uj << \fk]m by replacing e — > —muj 2 /k c g in 
equation (10). The asymptotic behavior is then given by 
log((5R(r))) ~ —r/l c {u) + iujr/v(uj), where 



I c (lj) = 



i cil 



mw|Im(fc e ^ 2 )| 



voir 



mRe(fc e 1 / 2 ) 



From Eq.Q one can prove the existence of two regimes, 
(i) For ui < uj* , the imaginary part 1/v(oj) = 0: one finds 
a pure exponential decay, with a characteristic length of 
order l c . (ii) For uj » uj* , the decay length decreases 
as I c {oj) ~ /(e) 1 / 4 ^ 1 / 2 , whereas the velocity of the cor- 
responding vibrations grows as viuj) ~ /(e) 1 / 4 ^ 1 / 2 . A 
similar behavior above uj* has been predicted by one of 
us for z > z c ^ 

We can also considering the case of our networks im- 
merse in a viscous fluid with viscosity 770 and without 
hydrodynamic interactions. The response in the over- 
damped regime is calculated by replacing the inertial 
term raw 2 by the viscous force —irjoui in equation ^ 
and then in the effective stiffness Q. In this limit 



the elastic moduli is proportional to k e g, whose real 
and imaginary parts gives the storage and loss moduli 
G' and G" respectively. Two different regimes are ob- 
served, (i) For uj « uj*' 2 m/r]Q, the storage moduli fol- 
lows G" - l 6 c uj 2 /f(e) and the loss moduli G" ~ 1\uj. The 
response is given by I c (uj) ~ l c and v(uj) ~ f(e)l~ 3 , dis- 
playing a similar exponential decay as the one observed 
in undamped system, but with propagating waves, (ii) 
For uj » uj* 2 m/rjo both, the storage and loss moduli 
are proportional to y f f(e)w, as also observed in Ref.^ 
The decaying length I c (uj) ~ /(e) 1 / 4 ^; -1 / 4 , while the ve- 
locity v(uj) ~ /(e) 1 / 4 ^ 3 / 4 . Our scaling for the charac- 
teristic length scale I c (uj) differs from the recent result 
of Tighe In our opinion, the difference stems from an 
incorrect definition of the length scale characterizing the 
elastic response in floppy materials (see footnote^ . 



4 Beyond mean field 
4.1 Fluctuations 

The obvious difference between the mean response to a lo- 
cal strain (Figj3j) and a typical one (Figj2]) indicates large 
fluctuations. A priori mean field models as the effective 
medium do not enable to capture those. However, it is 
possible to combine EM results with additional consider- 
ations to estimate the amplitude of these fluctuations. 

We denote the dipole of forces generated by 

changing the rest length of the spring ij by a distance 
one: 



131 



<j 



(12) 



where the stiffness coefficient k is set to unity. From the 
definition of the stiffness matrix, Eq.([2]), and Eq.(12), we 



t Tig hf m defines a length scale \f = f/F as the ratio between the 
typical contact force / carried by springs to the viscous forces F 
exerted by the fluid. F is proportional to the velocity of the parti- 
cles, and therefore to their displacements time the frequency u). It 
is thus readily extractable from our results. The contact forces / 
are such that forces are balanced on each node, which implies on 
average that the spatial derivative of the contact forces / are the 
external forces F. If one considers for example the response to a 
local perturbation, our result implies that (/}/|(F}| is of order of 
the length scale l c (t^) on which the mean elastic response decays. 
Thus this definition is consistent with our results. However a dif- 
ferent result is obtained if one considers the ratio W (f 2 ) / -J (F 2 ), 
as done by Tighe. These two quantities differ because the mean re- 
sponse is much smaller than the fluctuations around it (see section 
4.1 1. Therefore, associating the quantity = \J (f 2 ) / \J (F 2 ) to a 
length appears unjustified. In the zero frequency limit we can ac- 



tually calculate this last expression using the formalism developed 
inEH. For isotropic random floppy networks one gets (f 2 ) ~ u>* 
and (F 2 ) ~ oj'^ 1 , leading to a ratio ^/{f 2 )/{F 2 ) ~ w*" 1 ~ I 2 
consistent with the numerical results of Tighe. 



can write: 



M = ^2\F ij ){F 



(13) 



where the sum is taken over all the bonds. Note that any 
floppy mode |<5Ro) has by definition no restoring force 
.M|£Ro) = 0, implying that (5R |Fjj) = for all con- 
tacts ij. 

Thus, the response <5R^ to the elongation of a spring 
ij, which is equivalent to the response to a force dipole 
\Fij), has no components along floppy modes. Therefore 
the equation AiSRij = |P«) can be inverted. Using the 
spectral decomposition of A4, one gets: 




\SRij) = ^2\\SR u )(SR u \F iJ ), 



(14) 



cj>0 



where lu 2 and |5R W ) are the non-zero eigenvalues and 
the corresponding eigenvectors of the stiffness matrix 
M J 19 l 25 l Eq.(14) implies that the norm of the response 
follows (SRijlSRij) = E^oi^lF^} 2 /^. 
ducing the average amplitude of the response SR 2 



Intro- 



^2( ii \(S'Rij\S'Rij}, one gets using Eq.(13| 



SR 2 



Fig. 5: Rigid regions, shown as red diamonds and magenta 
square, are induced by freezing the nodes outside an external 
radius (red) and an internal radius (magenta) respectively. 
Circular nodes (blue) show the minimal floppy region. Freez- 
ing any extra particle rigidifies the entire system, (b) I*, as 
defined in the text, vs. Sz. Red line corresponds to P ~ l/5z. 



4.2 Pinning boundaries 
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For comparison, the total amplitude of the mean response 



can be calculated from Eq.(ll), and one obtains for the 
norm square J 5R m {r) 2 d d r ~ e^ _1 . Thus for d > 2, the 
norm of the mean response vanishes as 8z — > 0, whereas 
the norm of the fluctuations diverge. Thus relative fluc- 



tuations must diverge at small e, as observed in our data. 

The most simple scenario is that the fluctuations of 
the response (SR(r) 2 ) decays with the same characteris- 
tic length oc l c characterizing the mean response. Mak- 
ing this assumption, which is numerically verified (see 
below), and using the result of Eq.(15) , the angular av- 
erage must read: 



We now turn to the spatial properties of floppy modes. 
-((JR^I.MI^RcjThe response to the stretch of a spring, exemplified in 
Fig. |2]), is a floppy mode of the network where this spring 
is removed. We have thus shown that floppy modes can 
(15]be localized (in the sense of presenting an exponential de- 
cay) on a length scale l c ~ 1/VJz (with a log correction 
for d = 2). We now extend a previous counting argu- 
mentP^H to floppy networks (see alsc^Sl) and show that 
floppy modes are also characterized by another length 
scale I* ~ 1/Sz >> l c . Below I* strict localization is 
impossible: floppy modes cannot have a smaller compact 
support. Essentially, I* is the length at which the number 
of constraints r'j -1 that result from freezing the bound- 
ary of a system of size r/ is equal to the number of floppy 

d 



SR t (r) = (6R(r) 



rV2 



c (d-2 



» 4 h(r/l c 



(16) 



where log(h(x)) ~ — x for x » 1, and the e dependence 
is determined by Eq.(15). Rescaling the distance by l c 



leads to a very good collapse of the response at different 
coordination values, as shown in Fig. (|4j>d). Note that 
in two dimensions the amplitude does not need to be 
rescaled. 



modes in the bulk Sz . This counting argument suggests 
that if the boundaries are frozen on a scale I* or smaller, 
floppy modes must vanish. 

To test this prediction, we fix all the nodes outside a 
circle of radius r/. A frozen boundary induces a rigid 
region where floppy modes are forbidden. We determine 
this region using the pebble algorithm,^! as shown in 
Fig.|5]a). As rf decreases, all floppy modes eventually 
vanish. For any z we can define n* , the number of nodes 
involved in the last floppy mode, and I* = n* 1 / d = \pn? 
for d = 2. Our measurements are shown in Fig.([5]b) 
and follow the prediction I* ~ 1/Sz. This result implies 
that exponentially small displacement at distances r > 
l c cannot be neglected when the rigid-floppy transition 



induced by freezing boundary conditions is considered. 
On the other hand, our results imply that floppy networks 
that are stabilized by pinning boundaries on the scale I* 
present soft modes with exponentially small frequencies 
near the center of the sub-system, since there are floppy 
modes with displacements of very tiny amplitude near 
the boundaries, of order exp(— l*/l c ) ~ exp(— 1/y/Sz). 
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Fig. 6: a) Average displacement SR t vs. the distance from 
a local strain r for floppy networks plus weak springs (see 
text) with N = 90000 nodes and fc wcak /fc = 10~ 7 . b) Av- 
erage displacement SR t vs. the rescaled distance r/l c , using 
A(c p , c a ,r) = f.3. The response of strictly floppy networks 
(filled symbols) has been included for comparision. Inset: The 
initial decay magnified. 



5 Weak interactions 

As discussed in introduction, in some materials floppy 
modes are stabilized by weak interaction, as is the case 
in covalent glasses. In order to model these weak in- 
teractions, we consider (see for example^) floppy net- 
works with k — 1 and add a number of weak springs 
of stiffness fc wca k "C mLo* 2 , which gives a finite elasticity 
to the network. Floppy modes then gain finite frequen- 
cies, of order the characteristic frequency scale associ- 
ated with the weak interaction uj c = y/k^^Jm. A gap 
in the density of states remain present in the frequency 
range [ui c , u*], if the weak interaction is weak enough, 
ke. \/k W eak/ m << w*. In this limit the elastic mod- 
uli scale as /i wca k ~ k wca ^/8z^ where Sz is the excess 
coordination of the network of strong interaction. At 
long enough wavelengths the material, to good approxi- 
mation, must behave as a continuous elastic medium, and 
for cj << cj c the density of states must follow a Debyc 
behavior D{uS) ~ u) . Extracting a velocity of sound 
from the elastic moduli and computing the wavelength 
at w ~ w c , one obtains a wave length of order l c , which 



thus characterizes as well the length scale above which a 
continuum description becomes a good approximation. 

We now argue that concerning the response to a local 
strain, the limit fc wca k — ^ is singular: for any fc wea k > 0, 
the response decays as a power-law at large distances 
r >> l c , even in the limit fc WGa k — > 0. The singularity 
of this limit can be seen by decomposing the response 
to a local strain in two parts. First, we consider the 
displacement of characteristic amplitude SR that would 
follow such a strain in the absence of weak interactions, 
i.e. with fc wca k = 0. Second, we consider the additional 
displacement induced by the presence of weak springs. 
Indeed due to the weak springs forces are not balanced 
after step one, and forces of order F wca k ~ fc W cak<5-R have 
appeared on the nodes. In the limit fc wca k — > 0, these 
forces are vanishingly small and lead to no relaxation on 
the modes of non- vanishing frequency uj > lo* . How- 
ever the spectrum now presents modes (stemming from 
the floppy modes that exist as fc woa k = 0) of character- 
istic frequency uj c . These modes relax due to the un- 
balanced weak forces with amplitude x that must satisfy 
mu> 2 x = k wcak x ~ F woak ~ k wcak SR, implying that x is 
independent of fc wea k- 

Although localization is lost, the coordination contin- 
ues to influence the response to a local strain (see Figjf^i). 
In particular, near the imposed strain for r << l c , the re- 
sponse is weakly-influenced by the weak interaction, and 
decays with the same characteristic length, as appears in 

Fig(§5 . 

6 Rheology of dense suspensions 

6.1 Random networks under a global shear 

It is interesting to consider the behavior of the present 
networks under shear, if they were placed in a viscous sol- 
vent. This problem is formally related^ to the afflne sol- 
vent model of suspension flow of hard particles, where hy- 
drodynamic interactions between particles are neglected. 
We consider that our networks are made of point parti- 
cles, subjected to a viscous drag F proportional to the 
particles velocities V in the reference frame of the sol- 
vent, which is assumed to follow an afflne shear. One can 
show, from the definition of this model, that the viscos- 
ity r\ is proportional to the non-afflne velocity squared 
of the particlesJ 19 * 28 * 31 ! The viscosity is proportional to 
the ratio P/7 2 , where P is the total power dissipated 
and 7 is the strain rate. The power dissipated follows 
P ~ (F • V) ~ ||V|| 2 . It is convenient to write the veloc- 
ity field in terms of the non-afflne displacement following 
an infinitesimal strain SR/Sj, i.e. V = jSR/Sj. One 
find that r\ ~ ((8R/5j) 2 ) , which relates the viscosity to 



the non-affine response to shear. 

For isotropic floppy networks It was found numeri- 
cally^ that ||<5R/<5 7 || - 1/VSz, implying 77 - 1/Sz. Us- 
ing our previous result on the response to a dipole in 
floppy materials and a simple hypothesis it is straight- 
forward to derive this result, thus extending a previous 
derivation valid for z > z c j^ see also.^ Here we per- 
form this calculation for completeness, and because the 
present argument can be extended to predict that the cor- 
relation length under shear is l c . We consider an affine 
shear strain £7 applied on the network. After such a 
strain, unbalanced forces appear on the nodes: |F,5 7 ) ~ 
where the sum is taken over all the bonds 
and \Fij) correspond to a dipole of force as defined in 



( 12 ). The coefficients jij are equal to T-riij where T is 



the strain tensor, which is linear in Sj at first order. The 
displacement field can be written as a linear combination 
of responses to local dipoles |<5R,5 7 ) = Yluj) TijI^Rij), 
where |<5Rij) is defined in Eq.(14). Two-point displace- 
ment correlations in space obey 



C*«5 7 (r) = (5R 57 (x + r) • 6R Sy (x)) 

(ij)(ki) 



•r).5R w (x)). (17) 



We now make the assumption that in random isotropic 
networks (as those considered here), the response of dif- 
ferent dipoles is weakly-correlated (this assumption turns 
out to be incorrect for flow of particles where subtle cor- 
relations are present in the structures visited by the dy- 
namics). Using this assumption: 



C Sl {v) ~ <5 7 2 5^Rtf(x + r) ■ (JRtftx)), 



(18) 



(ij) 



where we used 7 2 ,- 
finds 



5j . Combining with Eq.(15) one 



C 57 (0) = (<5R 57 (x) • JRi 7 (x)) ~ 6~f 2 5R 2 d 



6j2 
e 



implying that the viscosity diverges as 77 ~ C<5 7 (0) /8j 2 



Moreover, ( 18 ) indicates that the correlation length in 



( 18 ) is essentially the length scale appearing in the re- 
sponse to a dipole l c . (JR^- (x + r) • JRy(x)) must vanish 
when r is larger than the length l c where the response to 
a dipole is localized. On the other hand, correlations are 
not expected to vanish on a scale much smaller than Z c , 
since the mean dipolar response will already give corre- 
lations on that scale: 



N ^ 

(ij) 



((SRy (x + r)) • ((5R y (x))dx ~ e~ 



6.2 ASM model of suspension flows 

We predict a plateau of modes above a frequency uj* ~ 
Sz. This plateau also appears in models of suspension 
flows.^ This observation supports that the plateau is 
generically present in floppy materials, and that the fre- 
quency response above uj* will be well described by the 
effective medium predictions derived here. There is, how- 
ever, a crucial difference between the spectra of our net- 
works and those of shear flows: in the latter one mode 
appears at a frequency w m i n « uj* that strongly couples 
to a global shear and leads to a divergence of the vis- 
cosity r\ ~ l/<5z 2 - 85 ,n23 which is sharper than in isotropic 
networks where rj ~ 1/Sz. Therefore, typical configura- 
tions visited in dense flows differ qualitatively from our 
random networks, for reasons explained in Ref.^ Never- 
theless, we expect some of our present predictions to hold 
in flow. 

For instance, the counting argument comparing surface 
effects and bulk degrees of freedom should hold, and I* is 
expected to characterize the effect of freezing boundaries. 

We do not know what the spatial correlations of the 
velocities should be in flows, because in this case the ve- 
locity is dominated by the lowest frequency mode that 
is not captured by the present analysis.! 19 * 32 * However the 
response to a local disturbance (such as the formation of 
a new contact) must couple predominantly to the modes 
above uj* , and is thus expected to affect the flow on a 
length l c . This can be shown using spectral decomposi- 
tion for the response under a dipole of force |F^). De- 
composing the spectrum in two parts (above and below 
uj*) Eq.dli) reads 



+ V ~|<5R w )<JR w |F y >. (19) 

^— ' UJ 

Numerical results show the appearance of a single mode 
below uj* however a theoretical estimation in the large 
TV limit supports that the distribution of modes between 
uj miD and uj* is given by AninM oc uj{uj 2 -uj 2 min )^ 2 ^ 2 ^ 
Following the argument of section |4.1| the average am- 
plitude of the response is given by 



— ^RylSRy)- ^ 
(ij) vj^l 



D min (uj) 



duj - 



D(uj) 



duj . 



I* 



The relative contribution follows I m i n /I* ~ Sz^ 1 ) 
(with a a logarithmic correction in two dimensions), in- 
dicating that in the limit Sz << 1 the contribution of 



the modes below u* becomes vanishingly small and can 
be neglected. The modes above w* in flow have the same 
density of states and are expected to have the same prop- 
erties than those of isotropic networks. Thus the response 
to a local disturbance in flow must decay exponentially 
with a typical length l c . This simple argument does not 
hold for the velocity correlation in flow: the modes in the 
plateau do not contribute significantly to the response to 
shear, which is rather dominated by the lowest frequency 
modes. 

7 Discussion and open questions 

We have argued that the elasticity of floppy networks is 
characterized by a gap in the vibrational spectrum, and 
by two length scales l c and I* that diverge near jamming 
(i.e. Sz — > 0). The existence of two lengths raises the 
question of which properties are governed by which scale. 
Our work supports that the length scale characterizing 
the response to imposed forces, and to most standard 
observations, is l c . On the other hand I* characterizes 
the effect of pinning boundaries. Such effects are subtle, 
and can depend on surprising ways on the type of elastic 
networks considered.^. 

For the strongly disordered floppy isotropic networks 
that we consider, one example of a question that remains 
to be explored is the evolution of elasticity (e.g. the shear 
modulus G) when boundaries are pinned at a distance L. 
For L > I* , the system remains floppy and G — 0. For 
L << l Cl one expects a mean-field argument to apply: 
pinning boundaries is equivalent to adding springs, lead- 
ing to an increase of coordination Az ~ 1 [L. For z > z c it 
is known that G ~ z — z c ,3Sand thus we expect G ~ 1/L. 
The behavior of G at intermediate length I* >> L >> l c 
remains to be explored. 

Finally we compare our results with previous works 
in amorphous solids made of soft repulsive particles, for 
which z > z c . Near the unjamming transition where pres- 
sure vanishes the coordination approaches the Maxwell 
threshold from above (z — > z c ). A vanishing frequency 
scale uj* ~ z — z c was predicted to characterize the 
low-frequency part of the spectrum,! 35 * 36 ! as confirmed 
numerically.^ The same theoretical argument^ indi- 
cated that boundaries affect elasticity on a length scale 
I* ~ l/(z — z c ). It was later argued,^ based on numer- 
ical observations, that I* characterizes the response to a 
point force in packings of particles. Another length scale 
l c was observed numerically to characterize the response 
at a frequency w*^and to affect transport,^ behaviors 
that were well-captured by effective medium.^ Our work 
extends these results to floppy materials with z < z c , 
where these lengths and frequency scales characterize the 



phonon gap and the localization of floppy modes. How- 
ever the comparison underlines an important disparity: 
for floppy networks we find both numerically and theo- 
retically that l c characterize the response to a local force, 
whereas I* appears to characterize the response to a point 
force in packings .^Sl More numerical and theoretical in- 
vestigations are needed to understand this difference. 
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